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First, we briefly review the conditions under which the entropy Sq can be extensive (Tsal- 
hs, GeU-Mann and Sato, Proc. NatL Acad. Sc. USA (2005), in press; cond-mat/0502274 ), 
as weU as the possible g-generalization of the central limit theorem (Moyano, Tsallis and 
GeU-Mann, cond- mat/0509229 ). Then, we address the q-triplet recently determined in the 
solar wind (Burlaga and Vinas, Physica A 356, 375 (2005)) and its possible relation with the 
space-dimension d and with the range of the interactions (characterized by a, the attractive 
potential energy being assumed to decay as r"" at long distances r). 

§1. On how can Sq be extensive 

Clausius introduced in 1865 the concept of entropy, S, without any reference to 
the microscopic world, i.e., to the possible microscopic configurations of the system. 
Indeed, at that time, most physicists did not believe in the existence of atoms and 
molecules. Boltzmann introduced, a few years later, the first connection of entropy 
to the microscopic world, namely 

w 



Sbg = -k'^Pilnpi . (1-1) 



i=l 

This expression is compatible with everything that Clausius expected for the entropy 
on thermodynamical grounds, in particular to be extensive. Indeed, if a system is 
composed of A'^ identical and distinguishable elements (or subsystems) that are prob- 
abilistically independent, i.e., '^j^^^^^tj^"~^^'^ = ''^i]^^^^ ■■■^tj^ ^ immediately verify 
that Sbg{^) = -^5'_bg(1)- Moreover, if the elements systems are quasi-independent 
(in a sense that needs of course more precise qualification, but which we do not 
address at this point), we have that Sbg still is extensive, i.e., SBGi^) oc N if 
N » 1. 

An interesting question arrives if the elements of the system are far from inde- 
pendence. Such is the case, for instance, when the elements have nontrivial (strictly 
or asymptotically) scale-invariant correlations, i.e., when the set of probabilities as- 
sociated with N elements is strongly correlated to the set of probabilities associated 
with (A^ — 1) elements. In such a case, which appears in fact quite frequently in 
natural and artificial systems, it can be shown that Sbg is not extensive. Instead, 
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the entropy 

^1 = ^ qit e 7^; Si = Sbg) ■ (1-2) 

can he extensive for a special value of the index q different from unity, as shown in.^) 
The most general binary case is given by {tt^^ , tt^^ } = {tt^^ , 1— 7rf ^ } if iV = 1, by 

Table I if A'' = 2, by Table II if A'' = 3, and so on. For arbitrary it can be thought 
as a A^-dimcnsional "hypcrcube" having W = 2^ probabilities {'^i^i^^'^'i^ '^'^^ } such 
that A1+A2+...+AM _ -I 
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Table I. Joint and marginal probabilities for N = 2, i.e., two binary subsystems Ai and ^42. 
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Table II. Joint probabilities T^t^^i^t^^'^^ for TV = 3. The quantities without (within) square-brackets 
[ ] correspond to state 1 (state 2) of subsystem ^3. 

In general, for A/" = 2, we have it^2^^^ 7^ ^^21^^^'^ foi" A/" = 3, we have 

41+A2+A3 / A1+A2+A3 / _Al+A2+A3 , A1+A2+A3 / A1+A2+A3 / „Ai+A2+A3. 

^112 T T^l2l T '•211 '^^^'^ ^122 T T^212 T ^221 ' 

and so on for A'' = 4,5,.... If we interpret A'' as the "time" of a time series, it 
is precisely these discrepancies which characterize "long-memory". We shall con- 
sider from now on the simple case where there is no "long-range" memory, i.e., 

^12 ~ ^21 = ^21 (witn TT]^^ = r2o and 1^22 = ^^22]) ^112 ~ 

A1+A2+A3 _ A1+A2+A3 ^ J A1+A2+A3 _ A1+A2+A3 _ /I1+A2+A3 ^ 

^121 — ^211 — ^31 ana 7r^22 — ^212 — ^221 — ^^32 

(with Ti^l^^'^^^''^ = and 71^2^"^^"^^^ = rss), and so on. This simple case can be 

represented on a triangle, as shown in Table III. 

We may impose on this triangle a very special correlation, namely the Leibnitz 

rule (see details in^^), i.e., rp^^n + TN,n+i = TN-i,n {n = 0, 1, A'' — 1; A'' = 2, 3, ...). 

This correlation makes the system scale-invariant in the sense that all the marginal 

probabilities associated with the A/'-system coincide with the joint probabilities of 

the {N — l)-system. The system may satisfy Leibnitz rule either strictly or only 
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asymptotically when N ^ oo. We will say respectively that the system is strictly 
or asymptotically scale-invariant. When the system is strictly scale-invariant it is 
enough to give one clement of each row (e.g., {rjsfo}, VAT) to fully determine it. If 
we choose rjvo = ^/{^ + 1) we obtain the Leibnitz triangle itsclf,^^ and can verify 
that it is "Boltzmannian" in the sense that it is only for q = 1 that we obtain 
extensivity. In other words, Sbg{N) (x N {N » 1) as can be seen in Fig. 1(a). 
We may choose a more complex set of probabilities that gives extensivity only for 
q = 1/2 as shown in Fig. 1(b). For brevity of space we skip here the details 
of this set, which is asymptotically scale- invariant. They details can be seen in.^) 
The main reason is, however, quite straightforward. The number of states whose 
probability differs from zero increases exponentially with N in the former case (Fig. 
1(a)), and only as N'^ in the latter (Fig. 1(b)) *\ It is then clear that the property 
SiAi+A2)/k = [Sq{Ai)/k]+[Sq{A2)/k] + [Sq{Ai)/k][Sq{A2)/k], which led to the term 
" nonextensive entropy" for Sg, is valid only under the (explicit or tacit) assumption 
of independence. 

Our overall conclusion at this point is that extensivity does not only depend on 
the specific functional form of the entropy but also on the composition law that we are 
using to form the total system out of its subsystems and their possible correlations. 
We have exhibited here systems of different nature (and different values of q) which, 
nevertheless, satisfy the Clausius requirement of being extensive. The various sets of 
probabilities that enter in the calculation of the entropy are, of course, expected to 
have basically a microscopic nonlinear dynamical origin. Strong chaos (i.e., at least 
one positive Lyapunov exponent for classical systems) would typically yield q = 1, 
whereas weak chaos (i.e., the maximal Lyapunov exponent vanishes) would typically 
yield q I, but the analysis of this fundamental point is out of the scope of the 
present short paper. 



(A^ = 0) 1 

{N = 1) no ni 

{N = 2) r2o r2i r22 

(N = 3) rso rsi r32 rsa 

[N = 4) r4o r4i r42 r^a ru 



Table III. Most general set of joint probabilities for A'^ equal and distinguishable binary subsystems 
for which only the number of states 1 and of states 2 matters, not their ordering. These proba- 
bilities satisfy J^^^q (N-!^n)'. n) ~ ^- particular simple case of independence corresponds 
to r^„ = (7r^)^-"(l-7r^)". 



*' If the number of states whose probability differs from zero increases as jj,^ with > 1 (as A'"'' 
with p > 0) the entropy Sq is extensive only for g = 1 (for q = 1 — 1/ p). 




§2. On the possible g-generalization of the central limit theorem 

In the case of the canonical ensemble we must optimize Sg under the following 
constraints^^ : 

w 



i=l 
W TP Ji 

= Ug, (2-2) 



J2Pi = h (2-1) 

i=l 

where Ei is the energy of the i*^ configuration and Uq some fixed quantity. The 
optimizing distribution of energies is given by p{Ei) = — — -/^^(g.-e/g) > Pq being an 
inverse temperature-like parameter and = [1 + (1 — q) (ef = e^). This ex- 

pression can be rewritten as p{Ei) = where (5q is a function of Uq. This 

corresponds to the case when the microscopic configurations are discrete. Its con- 
tinuous analogous corresponds, in the simplest case, to optimize Sq = /^LJl^^MflL 

with the constraints f dxp{x) = 1 and '^j f^^l^^jy,^^'' = Uq. A very elementary case 

corresponds to E{x) oc x^, hence p{x) oc e~^^^ , where is a positive constant. It 
happens that a considerable number of natural and artificial cases do exhibit this 
distribution. This would be neatly understandable if this form is the attractor, in 
the space of the distributions, under some quite generic type of composition of the 
random variables that are being averaged, for example if the variables have a spe- 
cial correlation. If this correlation is scale-invariant in a sense similar to the one 
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that we have discussed in the previous Section, the mathematical structure would be 
that of a g-generalized central limit theorem. The possibility of such a theorem has 
been addressed in several occasions in the literature/^ and was specifically analyzed 
in.*^) Strong numerical indications became recently available^) which g-generalize the 
celebrated de Moivre-Laplace theorem. Let us briefly review them. 

We consider a triangle such as that in Table III, and assume the Leibnitz rule to 
be strictly valid. To determine all the probabilities it is enough, as already mentioned, 
to provide, for instance, r^o. We shall adopt the g-product^^^ rjvo = 1/ [N2^~'^ — 
{N - 1)] 1/(1-9) (see details in^)). It has been verified that, for increasing N, the 
probability distribution approaches, after appropriate centralization, rescaling and 
symmetrization, a distribution p{x) oc e~^^ with K some positive constant and 

(?e = 2-i (g<l). (2-3) 

Let us introduce the definitions 

u{x) = 2 - X , (2-4) 
K^)^-, (2-5) 

X 

respectively corresponding to additive and multiplicative dualities (zy^(x) = fJ-'^ix) = 
X, \/x). We then see that Eq. (2.3) can be rewritten as Qe = i^lJ'{q)- The two trans- 
formations II and ly appear (isolated or combined) very frequently (e.g.,^^^) within 
nonextensive statistical mechanics. They can provide some hint for understanding 
the meaning of the g-triplet that was conjectured some time ago,^^^ and was recently 
confirmed by Burlaga and Vinasi'^-' with data on the solar wind sent by the spacecraft 
Voyager 1 from the distant heliosphere. 

§3. On the q-triplet 

The Q'-triplet is defined as {qsem Qreh Qstat) where sen, rel and stat stand for sen- 
sitivity (to the initial conditions), relaxation and stationary state respectively (see^^^^ 
for details) . Let us assume that we are addressing a classical many-body Hamiltonian 
system with two-body interactions given by an attractive potential energy which at 
the origin behaves smoothly and which, at long distances, behaves like (a > 0). 
It is reasonable to expect that the g-triplet depends on a and on the space-dimension 
d. It is in fact expected that it depends only on the ratio a/d.^^) 

Naturally, for a/d > 1 we expect Qsen = Irei = Qstat = 1- But for < a/c/ < 
1, the functional forms still arc unknown. We expect, however, qsenipt/d) < 1 < 
(lstat{ot/d) < qrei{c(/d). A simple possibility is as follows: 

u{a/d) , (3-1) 

fiia/d) , (3-2) 

(3-3) 



_ a _ 

qstat ~ 

d 

qrel = — = 

a 
a 

qsen ~ ~1 • 
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This set (see Fig. 2) satisfies the following relations 

qstat + Qsen = 2 , (3-4) 
Qsen Qrel — 1 ) (3'5) 
Qrel (2 - Qstat) = 1 • (3-6) 

Relations such as these might be applicable for Hamiltonian systems. 
Another simple possibility is as follows: 

ce 

qstat = 2 ~ = i^(«/<^) > (3-7) 

Qrel = - =ii{a/d), (3-8) 
a 

Qsen = 2 = un{a/d) . (3-9) 

a 

This set (see Fig. 2) satisfies the following relations 

Qrel + Qsen = , (3-10) 
Qrel {2 - Qstat) = 1 , (3-11) 
{2 - Qstat) {2 - Qsen) = I . (3-12) 

Relations such as these might be applicable for dissipative systems. 

It was determined in^^) that {Qsen, Qrel, Qstat) = (-0.6±0.2, 3.8±0.3, 1.75±0.06). 
Since important dissipative phenomena might be present in the solar wind, it seems 
reasonable to compare the experimental data with conjectures (3.7) to (3.12) (In 
any case, relation (3.3) is incompatible with the negative value observed for Qsen-)- 
If we do make such comparison, it appears as physically tempting to use the above 
possible expressions at their values corresponding to the Coulombian interaction 
within the solar wind plasma, i.e., a = 1 and d = 3, hence a/d = 1/3. An acceptably 
close theoretical set could then be {Qsen, Qrel, Qstat) = (—1)3,5/3). Although outside 
the observational error bars, we might not wish to exclude it at this early stage 
of this complex problem. If we wish, however, to get theoretical values closer to 
the observational ones, we might consider a set of relations such as (3.7-3.9) but 
somewhat more general in order to have a flexible scenario. We can then analyze 
the following more general forms: 

(3-13) 
(3-14) 
(3-15) 

where A, B and C are constants. We verify that, by construction, all g's equal unity 
for a/d = 1, as desired. For {A,B,C) = (2,0,2) we recover Eqs. (3.7) to (3.9). 

It is clear that a set {A,B, C) exists which reproduces the observational data. A 
different, more simple, possibility is to assume that the triplet is constituted by simple 
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Fig. 2. a/d-dependence of the g-triplet {qsen,<lrei,qstat)- Up left: set of Eqs. (3.7-3.9); Up right- 
set of Eqs. (3.16-3.18) Bo£tom;set of Eqs. (3.1-3.3). 



rationals. In that case a solution is provided by {qsem Qreh Qstat) = (~l/2, 4, 7/4) 
which satisfies^^ the following simple relations: Qstat = ^l^ilrel) cind qsen = IJ-^ilrel)- 
In this case we have {A, B, C) = (17/8, -1/2, 7/4), hence 



Qstat 
Qrel 
Qsen 



17 _ 9 a 

1 3d 

~2 + 2 a' 
7_ 3 d 

4 4 a ' 



(3-16) 
(3-17) 
(3-18) 
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Fig. 3. Solid and dashed lines correspond to Eqs. (3.7-3.9) and Eqs. (3.16-3.18) respectively. The 
NASA observational data^^' have been tentatively located at a/d = 1/3. 



This set (see Fig. 2) satisfies the fohowing relations 

Qrel "I" 2 Qsen 



1, 



{2qrel + l){n -Sqstat) 

27 

(7-4ggen)(17-8g^tot) 
27 



= 1, 



1. 



(3-19) 
(3-20) 
(3-21) 



The sets (3.7-3.9) and (3.16-3.18) are numerically not very different: see Fig. 
3. The former is mathematically extremely simple; the latter coincides with the 
experimental data within the error bars. We have not succeeded at this point to find 
further arguments for better understanding this complex phenomenon. Consistently, 
this entire Section cannot be considered as more than a very rough first discussion 
of the issue. 



§4. Conclusions 



We have illustrated that the extcnsivity of 5*^ for g 7^ 1 is possible in the presence 
of relevant correlations (e.g., strictly or asymptotically scale invariant ones) between 
the elements of the system. This seems to be consistent with the existence of a 
g-gcncralized central limit theorem, in the same way the standard, Gaussian, central 
limit theorem plays a crucial role in Boltzmann-Gibbs statistical mechanics. Some 
numerical indications already exist that point towards such theorem. These indica- 
tions provided a very simple relation based on additive and multiplicative dualities. 
This in turn enabled a first exploratory analysis of the values determined, for the 
g-triplet, from the data obtained for the solar wind at the distant heliosphere. 

One of the many interesting problems which remain open at this stage is the 
illustration of the g-generalized central limit theorem with gg > 1- A g-generalization 
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of the de Moivre-Laplace theorem that would lead to say qe = 3 — 2/q (q > 1) would 
certainly be very interesting. Indeed, such a relation implies that q increasing from 
unity to infinity makes qe to increase from unity to 3, the upper limit for gg-Gaussians 
to be normalizable. Also, q > 3/2 implies qe > 5/3, which corresponds to diverging 
variance of the ge-Gaussians. 
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